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Abstract 

The ATLAS detector at the LHC obtained the sample of p+Pb data at = 5.02 TeV with 

integrated luminosity of 25nb“^, which can be compared to the pp data obtained by interpolating 
pp measurements at ^/s = 2.76 TeV and 7TeV. Due to the excellent capabilities of the ATLAS 
detector, and its stable operation in heavy ion as well as proton-proton physics runs, the data allow 
measurements of the nuclear modification factor, ratio of heavy ion charged particle spectra divided 
by pp reference, in different centrality classes over a wide range of rapidity. The charged particle 
nuclear modification factor is found to vary significantly as a function of transverse momentum with 
a stronger dependence in more peripheral collisions. 
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1 Introduction 

The partonic structure of nuclei has been studied extensively in deep inelastic scattering experiments 
[1]. These measurements have established differences in the parton distribution functions (PDFs) in 
nuclei from those expected from an incoherent superposition of nucleons. These differences include the 
phenomena of shadowing, anti-shadowing [2] and the EMC effect [3]. These observations have inspired 
significant theoretical work (e.g. [4,5,6]), but to date the underlying physics of these phenomena is not 
completely understood. Proton-ion collisions at the LHC are capable of probing the nucleus at scales that 
have not been experimentally accessible before. This includes particle production from hard scattering 
processes with large transverse momentum transfer, q^, which may provide crucial insight into the large-g^ 
structure of the nucleus or even reveal previously unobserved phenomena. 

2 Analysis 

Nuclear modification effects are studied by comparing observables in heavy ion collisions with the cor¬ 
responding observables obtained in pp collisions. The rates of hard scattering processes in heavy ion 
collisions are enhanced relative to pp collisions due to the increased flux of partons per collision, i.e. the 
nucleus contains many nucleons capable of participating in a hard scattering. This enhancement depends 
on the thickness of the target nucleus, Tpb, and on the impact parameter of the projectile. It can be in¬ 
ferred from geometric models such as the Glauber model [7]. The particle yields in p-l-Pb collisions scaled 
by the average Tpb can then be compared directly to the cross section for the same process obtained in 
pp collisions: 


7?pPb(TT) 


1 1/A^evt d^Ap+pb/dydpT 

(Tpb) d^CTpp/dydpT 
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where RpPh is called the nuclear modification factor, 1/Aevt d^Ap+pb/dj/dpx is the per-event charge 
particle yield in p-|-Pb collisions measured differentially in transverse momentum (px) and rapidity (y), 
and d^cTpp/dydpx is the corresponding differential cross section measured in pp collisions. 

This analysis [8] uses p-|-Pb data collected with the ATLAS detector [9] in 2013 at ys))])' = 5.02 TeV 
with a longitudinal rapidity boost of 0.465 relative to the ATLAS laboratory frame. Further, the analysis 
uses pp data with i/s = 2.76 TeV recorded in 2013 and data with = 7 TeV recorded in 2010 and 2011. 
The p-|-Pb data have integrated luminosity (Lint) of 25nb“^, the pp data with ^/s = 2.76TeV and 7TeV 
have Lint ~ 4.0 pb“^ and 156 pb“^, respectively. 

The events used in this analysis are selected from several triggered samples. A minimum bias (MB) 
samples were obtained requiring signal in the Minimum Bias Trigger Scintilator. For jet triggered samples, 
jets were reconstructed in events that passed the MB requirement using the anti-kx algorithm [10] with 
the distance parameter R = 0.4. Events were selected by the jet trigger if they contained jets with Lx 
above a certain threshold. Multiple thresholds were defined. Only a fraction of all events which fired a 
trigger was randomly selected to be recorded for further analysis. This fraction was set differently for 



each trigger. Total spectra were obtained by merging spectra from MB and jet triggers, all of them scaled 
by reciprocal of the fraction to obtain MB spectra up to high px- 

The p+Pb event centrality classes are defined in terms of percentiles of the total transverse energy 
measured by the forward calorimeters in the pseudorapidity range 3.1 < 77 < 4.9 in the lead-going 
direction. The (Tpb) values are estimated from the Glauber model [11] for each of the studied centrality 
class. 

To study the effects of the detector response on the measurement, Monte Carlo (MC) simulation 
samples were produced using the PYTHIA event generator [12]. MC samples correspond to p-|-Pb and 
pp collisions of the same energy, the same detector configuration and the same beam direction. The MC 
samples were processed with the same algorithm as the data. For p-|-Pb, a sample of pp collisions at 
y/s = 5.02 TeV was simulated, but with the output boosted to match the boost due to the asymmetry 
of the collisions. This simulation output was overlaid onto minimum-bias p-|-Pb data to get events 
with generated jets and correct underlying events. Several samples of events were produced for different 
intervals of generator-level R = 0.4 jet px, besides the MB samples, in order to obtain good statistics 
precision over large range in charged particle px. The final pp (p-|-Pb) results were obtained by merging 
all pp (p-l-Pb) samples, while each sample was scaled by the cross section of the jet production in the 
appropriate jet px range. 

Charge particle tracks are reconstructed in the ATLAS Inner Detector [13]. Tracks are measured 
using a combination of silicon pixel detector (Pixel), silicon microstrip detector (SCT), and a straw tube 
transition radiation tracker (TRT), all immersed in a 2 T axial magnetic field. Charged particles typically 
traverse 3 layers of silicon pixel detectors, 4 layers of double sided microstrip detector, and 36 straws. 

The minimum px of tracks used in this analysis is 4 GeV. Each track is required to have at least 1 
hit in the Pixel detector, a hit in the innermost layer if such hit is expected by the tracking model, and 
at least 6 hits in the SCT. Additionally, tracks with px > 10 GeV are required to produce at least 8 hits 
in the TRT. These requirements select tracks with good px resolution and suppress the contribution of 
poorly reconstructed tracks. Yet they limit the analysis to the coverage of the TRT detector, thus all 
tracks are selected with j?;] < 2 . 

To ensure that the tracks originate from the event vertex, the transverse impact parameter, do, and 
Zosin0 (zq is the longitudinal impact parameter) are required to be less than 1.5 mm. Tracks selected for 
the analysis are also required to satisfy the conditions on the impact parameters significances: j do /cr^j, j < 3 
and j zo sin 0/azg sin e | < 3. The do and zq sin 9 parameters and their uncertainties are estimated by a vertex 
finding algorithm. 

To reduce the amount of fake tracks at high px, all tracks with px > 15 GeV are required to match 
anti-kx jets with the distance parameter R = 0.4. The main interest of this analysis are spectra of charged 
hadrons. Thus, electrons and muons coming from electroweak decays of heavy bosons are subtracted from 
the measured spectra. 

The raw charged particle spectra in p-|-Pb and pp collisions are corrected for fake tracks and secondary 
particles, for limited momentum resolution, and for tracking inefficiency. These corrections are functions 
of Px and rapidity in center-of-mass frame (y*) and are estimated using corresponding MC samples. 
First, the raw spectra are corrected for the fraction of tracks associated to secondary particles and for 
the fraction of tracks which cannot be associated to any particles in the MC, i.e. fake tracks. The 
fraction of secondary particles is at most 1% at px = 4 GeV. The fraction of fake tracks strongly depends 
on Px and at the highest measured px reaches 13%. It also depends on rapidity. Next, the resulting 
spectra are unfolded using the iterative Bayesian unfolding [14] to correct for the finite detector px 
resolution. The first iteration is sufficient, further iterations do not change the results by more than 2% 
even at the highest px. The spectra are then corrected for particle loss in the reconstruction by the track 
reconstruction efficiency. 

Once the differential pp cross section at ^/s = 2.76 and 7 TeV have been measured, the differential pp 
cross section at y/s = 5.02 TeV is obtained by interpolation. The interpolation is proportional to ln(v^) 
and it is performed for every px bin in each rapidity interval used in the p-|-Pb analysis. 

Systematic uncertainties are evaluated by varying individual sources within the ranges of their un¬ 
certainties and comparing the results to the results of the default analysis. Tracking and vertex pointing 
cuts contribute up to 5%. The uncertainties on the correction for fake tracks and on the correction for 
track reconstruction efficiency contribute up to 10% and 15% at the highest px, respectively. Uncertainty 
on the momentum resolution of the Inner Detector tracking system adds a significant contribution of 12% 
to the systematic uncertainties at high px. Difference up to 5% was found between parts of the run with 
opposite directions of the p and Pb beams in the LHC due to different jet trigger performance. The cal¬ 
culation of y* is made with the assumption that all particles are pions, this introduces an uncertainty of 


1%. The uncertainty related to the description of the inactive detector material is found to vary between 
1-7% depending on rj. 

The uncertainty on the event selection is accounted for in the uncertainty of the corresponding (Tpb) 
values. The luminosities of the pp data samples have an uncertainty of 3.2% and 4% for the measurements 
at ^/s = 2.76 and 7 TeV respectively. Possible distortion introduced by the interpolation algorithm reaches 
up to 12%. It is estimated as the difference between y/s and ln(ys) interpolations. For i?ppb, most of 
the systematic uncertainties between p+Pb and pp contribute to the final ratio, since the y* variable 
corresponds to different areas in the detector, as the two systems are in two different centre-of-mass 
frames. 


3 Results 

The final charged particle spectra are MB spectra (including scaled jet triggered spectra). The p+Pb yield 
is normalized per event, while the pp yields are normalized by corresponding luminosity and interpolated 
to the same center-of-mass energy as p+Pb sample. 

Left panel of Fig. 1 shows i?ppb in the rapidity interval —1 < y* < 1 for centrality class 0-90%. 
The figure shows unexpected increase of i?ppb with increasing px, compared to the lower px values. The 
i?pPb distribution reaches the maximum of approximately 1.4. This result is comparable with the CMS 
measurement [15] and both measurements show the same trend at high px. 

Right panel of Fig. 1 shows comparison of i?ppb with the interpolated pp reference and with a reference 
obtained by pure simulation at ^/s = 5.02 TeV. The increase of Rppb at high px is similar for these two 
distributions. 

Figure 2 shows the i?ppb distributions for 3 centrality classes (0-10%, 20-30% and 60-90%) and for 
the rapidity interval —2 < y* < 1.5. Especially in the 60-90% centrality class, there is a clear trend of 
increasing enhancement towards higher px. In the more central classes the enhancement is presented as 
well. 
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Figure 1 - The -RpPb distributions as a function of px in the 0-90% centrality class. Left: -Rppb measured for the \y*\ < 1. 
Right: -RpPb with interpolated pp reference (black points) and with simulated pp reference (red squares). Vertical bars 
indicate the statistical uncertainty of the measurement. Systematic uncertainties are plotted with step lines encompassing 
each point, the systematic uncertainties of (Tpb) are not included in the bands and shown with the gray-shaded bands at 
unity. 



































Figure 2 - The Rpph distributions as a function of px measured for the 60-90% (left), 20-30% (middle) and 0-10% (right) 
centrality class. Vertical bars indicate the statistical uncertainty of the measurement. Systematic uncertainties are plotted 
with step lines encompassing each point, the systematic uncertainties of (Tpb) are not included in the bands and shown 
with the gray-shaded bands at unity. 
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